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■ i-H . 1 INTRODUCTION 



ABSTRACT 

We construct a sample of 3,516 radio-loud host galaxies of active galactic nuclei (AGN) 
from the optical Sloan Digital Sky Survey (SDSS) and Faint Images of the Radio Sky at 
Twenty cm (FIRST). These have 1.4 GHz luminosities in the range 10 23 -10 25 WHz" 1 , 
span redshifts 0.02 < z < 0.18, are brighter than r* etro < 17. 77 mag and are con- 
strained to 'early-type' morphology in colour space (u* — r* > 2.22 mag). Optical 
emission line ratios (at > 3a) are used to remove type 1 AGN and star-forming galax- 
ies from the radio sample using BPT diagnostics. For comparison, we select a sample 
of 35,160 radio-quiet galaxies with the same r*-band magnitude-redshift distribution 
as the radio sample. We also create comparison radio and control samples derived by 
adding the NRAO VLA Sky Survey (NVSS) to quantify the effect of completeness on 
our results. 

We investigate the effective radii of the surface brightness profiles in the SDSS r 
and u bands in order to quantify any excess of blue colour in the inner region of radio 
galaxies. We define a ratio R = r e (r)/r e {u) and use maximum likelihood analysis to 
compare the average value of R and its intrinsic dispersion between both samples. R 
is larger for the radio-loud AGN sample as compared to its control counterpart, and 
we conclude that the two samples are not drawn from the same population at > 99% 
significance. Given that star formation proceeds over a longer time than radio activity, 
the difference suggests that a subset of galaxies has the predisposition to become radio 
loud. We discuss host galaxy features that cause the presence of a radio-loud AGN 
to increase the scale size of a galaxy in red relative to blue light, including excess 
central blue emission, point-like blue emission from the AGN itself, and/or diffuse red 
emission. We favour an explanation that arises from the stellar rather than the AGN 
light. 



Considerable uncertainties remain as to what controls the 
apparent link between the activity of active galactic nu- 
clei (AGN) during the time that a black hole is fed and 
star formation. The accretion disk surrounding a super- 
massive black hole (SMBH) emits highly energetic radia- 
tion and particles, and can for m powerful winds and/o r 
collimated, relativistic jets (e.g., Ivan Breugel et all [20041. 
The radiation and outflows might then affect the interstel- 
lar medium, triggering star formation which might be de- 
tectable as an excess of blue light from the central regions 
of galaxies. Alternatively, the activity of the AGN may be 
sparked by specific events in the galaxy's past. For exam- 
ple, there is morphological evidence that activity in radio 
galaxies mig ht be triggered by mergers an d galaxy inter- 
actions (e.g., iGonzalez Delgado et al.l 12006). which in turn 
could contribute to central blue light through enhanced star 
formation. In either case, we might expect an observable as- 
sociation between AGN activity, for which in this work we 
use radio loudness as a tracer, with bluer stellar continuum 
in the central regions of the galaxy. 

Pr evious work has hinted at such a relationship. For ex- 
ample, iMahabal et"al] j 19991 ) examined the central light of 



30 radio galaxies as compared to 30 normal galaxies from 
the Molongo Reference Catalogue and found excess blue 
light in the inner regions of the radio galaxies as compared 
with the control sample. This central excess is unrelated to 
the conv entional colour gradien t of th e broader distribution 
of stars. iGonzalez-Perez et al.l |2010l 'l studied colour vari- 
ations within galaxies fr om the Sloan Digital Sky Survey 
(SDSS; FYork et all l200ch DR7, finding marginally steeper 
colour gradients in massive galaxies with nuclear activity 
and concluded that this is due to a higher fraction of young 
stars in their central regions. 

In the p r esent paper we expand the work of 
IMahabal etafl |l999h . by combining optical information 
from SDSS with t he FIRST and NV SS radio surveys 
jBecker et al.lll995l ;l Condon et al.l fl998) to construct a lo- 
cal (0.02 < z < 0.18) sample of powerful radio- loud AGN 
host galaxies (R-AGN), and a control sample of 'normal' 
early-type ellipticals. We define R as the ratio of SDSS r* 
to u* band de Vaucouleurs effective radii and compare its 
value between the samples. 

The paper is structured as follows. In Section [2] we 
briefly describe the surveys used to construct our samples. 
Separation into the radio-loud and control samples is dis- 
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cussed in Section [3] Section [4] outlines the maximum likeli- 
hood analysis used to compare the average values of R for 
the two samples. In Section [5] we present a discussion of our 
findings in the context of radio-AGN fuelling mechanisms 
and the properties of the host galaxies. Throughout, we use 
a fiat ACDM cosmology with Q mo = 0.3 and Oao = 0.7. We 
adopt H = 70 kms" 1 Mpc -1 . 



2 DATA SOURCES 
2.1 SDSS 

SDSS is an imaging and spectroscopic survey covering 
~10,000 deg 2 , primarily in t he northern hemispher e. Its 
7th data release (SDSS-DR7; lAbazaiian et al.l 120091 ) con- 
tains over 350 million entries, of which about one million 
are confirmed galaxies with spectroscopic follow up. 

In the photometric sample, flux densities are measured 
simultaneously in five broadband filters (u,g,r,i, z), with 
effective wavelengths 3551, 4686, 6165, 7481 and 893lA 
(Fuku gita et al.l I1996T). All mag nitudes are given on the 
AB„ system (|Oke fc Gunnll 19831), and are determined using 
Petro sian apertures (see lBlanton et al.ll200ll ; IGraham et all 
120051 . for a comple te definitio n of th e Petrosian system) . 
In common with llvezic et al.l (|2002l ). we refer to SDSS 
measured magnitudes as u*,g*,r*,i* and z* due to the 
uncertainty in the absolute calibration of the SDSS pho- 
tometric system, which ha s been assessed at <0.03 mag 
(|Stoughton fc LuptonI 120021') . For a complete guide to the 
photo metri c system. see | Fukugita et al.l (| 19961 ). Gunn et"al] 
|200ll ) and lStoughton fc LuptonI |2002i ). 



2.1.1 Main galaxy sample (MGS) 

We use the main ga laxy survey spectroscopic sample (MGS, 
IStrauss et al. 2002). The MGS algorithm selected extended 
sources brighter than r* etro = 17.77 mag with Petrosian 
half-light surface brightness / X50 ^ 24.5 mag arcse c -2 , pro- 
viding ~ 90 galaxies deg -2 (Bl anton et al.l [200lT ) . Candi- 
date galaxies were then observed spectroscopically. 99.9% 
of the resulting redshifts have velocity errors < 30 kms - . 
The survey is unbiased, except for galaxies with close com- 
panions. Around 6% of galaxies satisfying the photomet- 
ric target criteria were not included for spectroscopic follow 
up due to a companion galaxy within the 55" minimum fi- 
bre separation, although some of these galaxies have been 
subsequently observed. The MGS contains ~680,000 spec- 
troscopically cojrfi£med_gala2des. The full target selection is 
described in lStrauss et al.l ([2002]). 

We corrected for Galactic extinction using the SDSS 
'reddening' parameter, which was derived from maps of 
the infrare d emission from dust across the s ky in accor- 
dance with lSchlegel. Finkbeiner fc Davis] (1 19981 ). We applied 
fc-correction s using the S PSS-derived 'kcorr' parameter as 
detailed in iBlanton fc Roweisl |2007l ). Typical corrections 
k u , k g , k r were (0 - 0.4), (0 - 0.3), and (0 - 0.1) mag. 



2.2 FIRST 

The FIRST survey l|Becker et al.ll 19951) utilized the VLA (at 
1.4 GHz in B array) to map the radio sky over ~9000 deg 2 in 



the Northern hemisphere and in a 2°. 5 wide strip along the 
celestial equator. It has ~ 8400 deg 2 overlap with SDSS, con- 
tains ~ 97 sources deg -2 at the 1 mjy survey threshold and 
reaches an rms sensitivity of ~0.15 mjy beam -1 . In this con- 
figuration, the VLA has a synthesized beam of 5.4" FWHM, 
providing accurate flux densities for small-scale radio struc- 
tures, but underestimating the flux densities of sources ex- 
tended to several arcminutes. At the 1 mjy survey threshold, 
individual sources have 90% confidence astrometric errors 
< 1". We adopted a mean spectral index of a = 0.7 (where 
Su oc ^ -a ) and obtained rest-frame 1.4 GHz power densities 
by applying k-corrections assuming the usual form 



L u , rext = (l + 2) (a-1) S„47r.Di 
where Dl is the luminosity distance. 

2.3 NVSS 



(1) 



The NRAO VLA Sky Survey (NVSS. ICondon et al.l [l998h 
mapped the radio sky (at 1.4 GHz in D array) north of —40° 
declination. The survey is complete down to a point source 
flux density of ~2.5mJy. The synthesized beam of 45" is 
much larger than FIRST, providing more accurate flux mea- 
surements for highly extended sources. Astrometric accuracy 
ranges from 1" for bright sources to 7" for the faintest de- 
tections. The entire survey contains over 1.8 million unique 
sources brighter than 2.5 mjy. fc-corrections are applied as 



3 SAMPLE SELECTION 

3.1 Cross-correlation: FIRST-SDSS 

We initially derived a FIRST-SDSS sample of radio galaxies 
and a comparison control sample with no nearby detectable 
radio counterpart. Figure[T]shows the distribution of angular 
separations between SDSS objects and their nearest FIRST 
counterpart. We find that at ~ 5" the distribution becomes 
dominated by random matches. Hence we adopt a match cri- 
terion of 2" angular separation to avoid such false matches. 
Our cross-matched catalogue contains 25, 931 galaxies and 
we denote this as our preliminary 'radio' sample. All un- 
matched sources are initial candidates for our control sample 
of radio-quiet galaxies. 

We define our efficiency as the fraction of matches in 
our radio sample which are physically real (including con- 
tamination from line of sight false matches, which cannot 
be accounted for here). In order to evaluate the number 
of random false matches in our sample, we offset the RA 
and DEC of the FIRST sources by 1° and re-matched to 
our set of 680, 056 SDSS galaxies, selecting objects within a 
2" radius. We found on average 58 random matches, corre- 
sponding to 0.3% contamination of the FIRST-SDSS match 
sample (>99% efficiency, Figure [T] red dotted line). 

We also need to know the completeness of our match- 
ing criteria - the fraction of correct matches we recover. 
A factor affecting this is th e possible exclusion o f large 
lobe-dominated radio sources l|Kimball fc Ivezidl2008l ). Both 
lobes will be included in FIRST but may be excluded in the 
cross-matched subset for lobe-core distances > 2", if the core 



The host galaxies of radio -loud AGN: colour structure 3 



2500 



2000 



1500 



- 000 



500 




.0 -0.5 0.0 0.5 
Separation, log,„(arcsec) 



1 .0 



Figure 1. The distribution of angular separations of the matches 
between FIRST and SDSS, truncated at a 20" radius (black). The 
distribution of matches between the FIRST positions offset by 1° 
and SDSS is shown as a red dotted line. This indicates the number 
of likely false matches. The net distribution of 'real' matches is 
shown by the blue dashed line, and our 2" cut indicates 99.9% 
efficiency of the matched sample. 



is weak. We estimated ~ 8% of real matches between the 
SDSS subset and FIRST that are excluded from our radio 
sample, by assuming the radio lobes are outside a 2" radius 
from the optical core, but within a 5" radius (see Figure 
[T}. Our m atching algorithm e fficien cy and completeness are 
similar to iKimball fc Ivezid l|200ct l. who match FIRST to 
SDSS-DR6 positions within a 2" radius for sources with z ^ 
2 (95% efficiency and 98% completenesf]) . 

However, we could not account for sources with FIRST 
lobes > 5" from the optical core or the population of 
FIRST double-lob ed sources with no detectable radio core. 
Ilvezic et all (|20Q2h cross-correlated all FIRST sources with 
SDSS, then identified potential double-lobed radio sources 
with undetected coreqj, estim ating these contr ibute less 
than 10% of all radio sources <|Best et all l2005al . estimate 
~ 5% of their radio-loud AGN sample had no FIRST detec- 
tion). 

Becker et all (|l995h estimated the FIRST catalogue to 
be 95% complete at 2mJy and 80% complete down to the 
survey limit of 1 mjy. Figure [2] shows integrated radio power 
against redshift for the preliminary radio sample (25,931 
sources) . The solid black lines trace the 1 mjy and 2 mjy 
peak flux density thresholds, the dotted and dashed lines 
show cuts we applied in redshift and radio power (see H3.3I 
and i]3.4p . and the red points are the resultant selection of 
radio sources. It is noted that for sources that are not well- 
described by an elliptical Gaussian model, the integrated 
flux density as derived by FIRST may be an inaccurate mea- 
sure of the true value. Such sources with radio powers cor- 



1 IKimball fc Ivezid (2008) estimate their completeness and effi- 
ciency by fitting a gaussian to the nearest-neighbour distribution 
(representing physical matches) plus a rising linear function (rep- 
resenting random matchings). 

2 Via comparing the mid-points of FIRST pairs to SDSS sources 
within a separation < 90" and accepted all matches with offsets 
< 3". 
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Figure 2. Integrated FIRST radio power vs redshift for the pre- 
liminary 'radio' sample (see text). The solid black lines trace peak 
radio power at the 1 mjy detection threshold and at 2 mjy, the 
dashed lines denote the redshift limits imposed in i|3.3l and the 
dotted line denotes the radio power cut imposed in i|3,4l The red 
points are the remaining radio sample after these cuts are applied. 



responding to flux densities below the 1 mjy threshold can 
be seen in Figure [5] 

We note that although the two cuts improve the sample 
completeness (all red points are above lmjy), many of our 
sources are below the 2mJy line, and so our sample com- 
pleteness cannot be > 80%. Combining this with our es- 
timation of completeness in our selection criteria (~ 90%), 
we estimate the actual completeness of our preliminary radio 
sample to be ~ 72% for sources brighter than 1 mjy. The in- 
completeness is dominated by statistical effects rather than 
the 10% effect of sources missing due to weak cores. 

To produce a preliminary control sample, we select 
sources from the SDSS subset which fall within the same 
spatial region as FIRST (~ 8000 deg 2 ). We selected all ob- 
jects from this subset which are not matched to a FIRST 
source within 2" of their optical core, providing ~625,664 
radio-quiet galaxies. The efficiency of this sample, based on 
Figure [T] is > 99.9% (i.e., < 1% contamination by radio 
sources detectable at the FIRST flux density limit) and the 
sample completeness is > 92%. 



3.2 Cross-correlation: FIRST-NVSS-SDSS 

For FIRST sources substantially larger than 5" some flux 
density is resolved out, leading to an increase in the survey 
threshold for extended objects and flux density underesti- 
mates for larger objects (|Lu et al.ll2007l ). Our FIRST-SDSS 
radio sample is large due to the number of FIRST candi- 
date sources near the 1 mjy survey limit. Its limitations are 
potential flux density underestimates, and incompleteness 
to extended radio sources, notably extended, double-lobed 
sources with no detectable radio core within several arcsec- 
onds of its optical counterpart. 

To attempt to quantify these limitations, we created a 
complementary sample of r adio-loud AGN ho st galaxies us- 
ing both NVSS and FIRST. llBest et all |2005al ) note that the 
45" resolution of NVSS is large enough that ~ 99% of all 
radio sources are contained within a single component, al- 
lowing for a higher sample completeness. However, NVSS is 
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Matched 
surveys 


Slim(mJy) 


Objects 


Completeness Efficiency 


FIRST- 


1.0 


25,931 


72% > 99% 


SDSS 








FIRST- 


2.5 


5,719 


99% > 95% 


NVSS-SDSS 









Table 1. Radio flux density limits, Sii m , of both the radio sam- 
ple (FIRST-SDSS) and the candidate CR sample (FIRST-NVSS- 
SDSS), the number of objects in each and the estimated com- 
pleteness and efficiency. 



less deep than FIRST and consequently the matched sample 
contains fewer galaxies. The NVSS-FIRST-SDSS matched 
sample therefore does not benefit from such good statistics. 

We used the Unifie d Radio Catalogue constructed by 
iKimball fc Ivezid l|200Sh (see sample 'C in their table 8), 
selecting sources which are detected by both FIRST and 
NVSS, matched to within 25". This selection yields a radio 
flux density limited (Si.4GHz > 2.5 mjy) catalogue (here- 
after FIRST-NVSS) containing 141,881 sources. We cross- 
matched these to the ~ 680, 056 MGS sources adopting the 
same 2" match criterion used in H3.ll Our cross-matched 
catalogue contains 5,719 galaxies and we denote this as the 
preliminary 'comparison radio' (he reafter CR) sampl e . Inte - 
grated flux densities as derived bv lKimball fe Ivezid (|2008l ) 
are adopted in our analysis. 

We estimated our matching criterion of 2" to be > 99% 
complete, and > 99% efficient (8 random matches were 
found when the FIRST-NVSS RA and DEC were offset by 
1°). The FIRST-NVSS catalogue used to create the CR sam- 
ple is 99% complete and matched with 96% efficiency (see 
table 2 of IKimball fc Ivezid [2008). We therefore estimated 
the comparison radio sample to be > 99% complete and 
> 95% efficient (see Table[T]). 5 of the 5,719 galaxies in the 
CR sample reside in the control sample derived in H3.ll and 
these were removed. 4,935 of the CR sample (86%) are also 
in the FIRST-SDSS sample. 



3.3 Redshift range 



In the following sections f i)3.3l3.8|) we discuss the secondary 
selection criteria applied to the FIRST-SDSS sample. The 
CR sample follows the same path, and a summary of the 
final CR sample properties are presented alongside those of 
the final radio sample in H3.9I 

We imposed a redshift cut of 0.02 < z < 0.18 to ensure 
that the galaxy spatial structure could be well examined. 
For z < 0.02, redshift is not a reliable distance indicator. 
Figure [3] shows less than 2% of matched sources have red- 
shift < 0.02. The MGS magnitude limit of r* etro < 17.77 
corresponds to L r * > 2.6 x 10 23 WHz -1 at z = 0.18, and 
our sample should be complete to this optical luminosity 
density. 

We restricted the SDSS-derived redshift confidence 
(zconf idence) to be greater than 95%, which cuts out 
12% of objects from the preliminary radio sample and 10% 
of objects from the preliminary control sample. After red- 
shift selection, 499, 418 radio-quiet galaxies remained in the 
control sample, ~70% of the parent MGS, whilst the ra- 
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Figure 3. Distributions of FIRST-SDSS radio sample redshifts 
for log 10 Pi.4GHz luminosity bins. We cut at 0.02 < z < 0.18, pro- 
viding a high sample completeness (> 98%) down to Pi.4GHz ~ 
10 23 WHz -1 . 



dio sample contained 17, 024 galaxies. The average error in 
redshift for both samples is < 0.002 and the redshift distri- 
butions of these two samples were indistinguishable at the 
1% level on the basis of a two-tailed K-S test. 



3.4 Radio Power 

iBest et all (|2005ah suggested that typical radio-lo ud AGN 
have powers above 10 23 WHz _1 at 1.4 GHz and ICondonl 
(1989) shows that Pi. 4 ghz = 10 23 WHz _1 separates the 
spiral starburst population from the AGN-E/S0 popula- 
tion in local galaxy fields. We therefore restrict the radio 
sample to galaxies harbouring radio sources with Pi.4GHz > 
10 23 WHz -1 , which selects 5,119 galaxies (30%) from the ra- 
dio sample. Figure[3]shows that this threshold in radio power 
is well matched to our redshift selection. We name this set 
of 5,119 galaxies the 'radio-loud' sample. The mean redshift 
of the radio sample increased from ~ 0.095 to ~ 0.136 when 
the radio power cut was applied. 

AGN can be split on the basis of large-scale radio struc- 
ture into FRI type, where radio brightness decreases out- 
wards from the cent re and FRII type, with edge-brightened 
lobes and hot spots (|Fanaroff fc Rilevll 19741 ) . FRIIs are more 
luminous Pi 78 mhz > 1.3xl0 26 WHz -1 and are rare in a low- 
redshift sample such as ours (Figure [3]), since this 178 MHz 
power corresponds to Pi.4GHz ~ 3.1 x 10 25 WHz _1 for a 
typical spectral index of a — 0.7. 



3.5 Removing type 1 AGN 

In unified AGN models (e.g. I Antonuccil 19931 ) the appearance 
of the central black hole and associated continuum of an 
AGN differ only in the viewing angle at which it is observed. 
Sources viewed face on (type 1) show broad emission lines 
that are absent in those observed edge on and where the 
broad emission line region is obscured by a dusty torus (type 
2) . Type 1 AGN are excluded from this study, as the optical 
continuum can be dominated by relativistically boosted non- 
thermal emission, which may overwhelm measurements of 
the host galaxy's properties. 
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The SDSS spectral classification pipeline automatically 
flags and excludes quasars from the MGS, but we chose to 
verify its reliability and the relative numbers of type 1 AGN 
re maining in our samp le. We followed the method outlined 
by iMasci et~all (|2010r i to identify broad line emission. Ha 
or H/3 emission lines exceeding 1000 km s^ 1 (FWHM) with 
a S/N > 3 and Ha/H/3 EW > 5 A were classified as broad 
line. Galaxies with both broad Ha and H/3 emission were 
classified as type 1 AGN. 3 galaxies of the 5,119 radio-loud 
sample have broad Ha and H/3 emission lines and were re- 
moved from the sample. 

184 galaxies (4%) of the remaining 5,116 radio-loud 
sample have broad Ha emission, but do not have broad H/3 
emi ssion. 

lOsterbrockl (|l989h classified these objects as type '1.9' 
AGN, which have substantial but not complete obsc u ration 
of the central continuum source. iKauffmann et al.l (2003) 
determine that the contribution to the observed continuum 
is not significant in these sources, so we retain these within 
our sample. The control sample contained 10 type 1 AGN 
and 1835 type '1.9' AGN. We removed the type 1 AGN to 
leave 499,408 galaxies within the control sample. 



3.6 Removing star-forming galaxies 

Type 2 AGN have narrow permitted and forbidden lines 
and their stellar continuum i s ofte n similar to normal star- 
forming galaxies. lYun et alj |2001) show a tight correlation 
between far infra-red luminosity (indicative of star forma- 
tion) and Pi. 4 GHz- This will cause a level of contamination 
by star- forming galaxies if Pi.4GHz is used as the sole tracer 
of galaxies hosting a radio- loud AGN. We should therefore 
remove the small subset of radio-loud galaxies in which the 
radio power arises from star formation (SF) and not from 
an AGN. 

AGN separation from SF galaxies in the local Uni- 
verse ca n be achieved via optical emission line ratio diag- 
nostics l|Baldwin. Phillips fc Terlevichl [l98ll . herein BPT). 
Emission-line ratios probe the ionizing source: for AGN, 
non-thermal continuum from the accretion disc around a 
black hole and in star-forming galaxies (SFGs), photoion - 
ization via hot massive stars. However, iBest et all {2005b) 
find no correlation between a galaxy being radio-loud and 
whether it is optically classified as an AGN. In agreement 
with this, we found no correlatiorff] between radio flux at 
1.4 GHz and the optical line ratios [Nil] /Ha or [OIII]/H/3 in 
our sample. Hence, a substantial fraction of radio- loud AGN 
would not be selected using BPT diagnostics, and were we to 
apply them to the radio sample it would be biased towards 
radio galaxies with particularly strong optical emission lines. 

We instead remove galaxies which are strongly iden- 
tified as non-AGN, i.e. star-forming. Despite this method 
leaving a small fraction of star-forming galaxies which are 
faint in optical line emission, all radio-loud AGN, whether 
optically bright or otherwise, will remain in the sample. This 
decrease in efficiency of the sample is preferential to a dras- 
tic decrease in completeness. A similar problem was identi- 



3 Spearman's rank correlation result is ±0.06 or less between the 
radio-flux and optical emission-line flux ratios for our radio galaxy 
sample. 
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Figure 4. BPT diagram for galaxies within the radio-loud sample 
with emission line ratios. 1,069 galaxies from the 5,116 radio-loud 
sample have emission line ratios and are plotted as grey crosses 
where the size of the symbol does not indicate the sizes of the 
error bars, which vary wi dely. The dashed curve ( eq . [2l indicates 
the demarcation given by [Kauffman n et al ] (|2003h between opti- 
cal AGN (above the line) and SF galaxies (below the line). 766 
galaxies plotted have emission line ratios with > 3<r significance. 
The density of these galaxies is shown by contours, 85% of the 
766 galaxies lie above the line as optically classified AGN. The 
remaining 15% (118 galaxies) lie below the curve (red circles) and 
we classified these as star-forming. We also show the maximum 
uncertainty for a point on the Kauffmann separator with S/N of 
3 (dotted lines). 69 objects detected at S/N> 3 (black crosses) 
lie between the Kauffmann separator and the upper maximum 
uncertainty line. 



tied bv lSadler etal] (|2002f ). who defined a sample of radio- 
loud AGN from the 2dFGRS catalogue. They found approx- 
imately half of the sample have absorption spectra similar 
to those of inactive giant ellipticals, and therefore would be 
mostly missed by optical AGN emission-line selection. 

In classifying galaxies as AG N or star-forming, weuti- 
lized the demarcation criterion of [Kauffmann et al.l ((2003) 

log([OIII]/H/3) > 0.61/{log([NII]/Ha) - 0.05} + 1.3 (2) 

plotted as the dashed line on Figure [4] 

Figure [4] shows the standard line ratio diagnostics for 
galaxies from the radio-loud sample with all four emission 
lines catalogued in the SDSS (grey crosses, 1,069 objects 
of the 5,116). 766 of these 1,069 galaxies have both optical 
emission-line ratios at S/N ^ 3 and their density on the 
BPT diagram is shown as contours. 118 (11%) lie below the 
demarcation line and are marked with red circles (112 out 
of these 118 have all four emission- lines with S/N ^ 5). We 
removed these 118 optically selected SF galaxies from the 
radio-loud sample, to leave a sample of 4,998 predominantly 
radio-loud AGN hosts. 

We then estimated the residual contamination expected 
in the radio AGN sample from star-forming galaxies by plot- 
ting the maximum uncertainty for a point on the Kauffmann 
separator with S/N of 3 (FiglU dotted lines). 69 objects de- 
tected at S/N > 3 (black crosses) lie between the Kauffmann 
separator and the upper maximum uncertainty line. There- 
fore the contamination expected in the radio AGN sample 
from star- forming galaxies is ~ 11%. 
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Our radio sample contains predominantly FRI galax- 
ies, which are usually hosted by giant elliptical galaxies and 
on average have weak or no optical nuclear emission lines 
|Lin et al.l |2010| . and references therein). Within our sam- 
ple of radio-loud AGN, 79% of objects do not have optical 
emission line fluxes. We estimate m 11% contamination of 
non-AGN (e.g. SFG) if galaxies without SDSS emission lines 
are similar to galaxies with bright line emission. 5 galaxies 
without emission line fluxes > 3cr lie below the demarca- 
tion line in Figure [4] and are potentially star forming, but 
without reliable line information we retain these within our 
sample. 

As discussed bv lBest et alj (|2005bf ). a potential short- 
fall of spectral classification of emission-line radio-loud AGN 
is that emission-line AGN activity is often a ccompanied by 
star formation fe.g.. [Kauffmann et al.| [2003y This star for- 
mation will give rise to radio emission, even if the AGN 
itself is radio quiet. For these sources, the optical spectrum 
could still be dominated by a (radio-quiet) AGN leading to 
clas sification as a n emi ssion-line radio- loud AGN. 

iMoric et ail \201(j ) (hereafter M10) derive a popula- 
tion of sources mat ched from all NVSS ga laxies in the Uni- 
fied Radio Catalog (|Kimball fc Ivezidl200g| ) . the SDSS-MGS 
and IRAS data. The matched sources are divided into star- 
forming, composite and AGN using standard BPT diagnos- 
tics. Star formation rates are derived via broad-band spec- 
tral fitting to the NUV-N1R SDSS photometry, and the aver- 
age fractional star formation/ AGN contribution to the radio 
power is estimated (see their table 2). They find that in 203 
composite galaxies, 81.3% of the total radio power is due 
to star formation. The variation of this fraction with radio 
power is not specified. 

Following M10, we defined 206 composite galaxies in 
our own radio sample (27% of the emission-lin e galaxies con- 
firmed at S/N ^ 3), using the dia gnostics of iKewlev et al.l 
(|200lh and lKauffmann et all (2003). If the average fractional 
contribution is independent of total radio luminosity, then 
we estimated an upper limit of ~ 27% of our radio sam- 
ple may have radio power boosted by star formation but 
possess a radio-quiet AGN. However, < 14% of the compos- 
ite sample defined by M oric have log[Pi.4(W Hz -1 )] > 23. 
iMauch fc Sadler! (|2007h find SF galaxies tend to have me- 
dian log[Pi. 4 (WHz~ 1 )] = 22.13, whereas AGN have a me- 
dian log[Pi. 4 (WHz- 1 )] = 23.04. Therefore, our luminosity 
cut will have significantly reduced the numbers of galaxies 
where star-formation is the principal contributer to the total 
radio power, and we expect far less than 27% contamination 
from radio-quiet, optically-loud AGN. 

We also cannot account for the population of 'compos- 
ite' radio-loud host galaxies with ongoing star formation 
that have been lost from our sam ple through this selection 
technique. IMauch fc Sadler! (|2007T ) estimate that ~ 10% of 
local radio sources may have a starforming spectrum but 
have radio flux densities dominated by a radio-loud AGN. 

3.7 Colour bimodality 

As low-redshift radio-loud AGN are hosted predominantly 
by elliptical galaxies, we attempted to constrain the control 
sample to contain only early-type galaxies. The colour dis- 
tributions of galaxies have been shown to b e highly bimodal 
jYan et al.ll20ol.l2010l ; lstrateva et alj|200 ll) . Figure [5] shows 
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Figure 5. Colour-colour plot for the entire control sample (grey) 
and 34606 visually classified ellipticals from Galaxy Zoo-control 
sample match (logip contours). Th e dashed line is u* — r* ^ 
2.22 mag from lStrateva et al. I ll200ll) and defines a morphological 
colour separator. 80% of the Galaxy Zoo-MGS sample lie above 
the line. The u* , g* and r* magnitudes are derived from SDSS 
columns: petrocounts — reddening — kcorr. 

a colour- colour plot in g*-r* against u*-g* as explored by 
IStrateva et aD (|200lh who define an optimal colour sepa- 
ratoiQ between early and late types of it* — r* ^ 2.22 mag. 
Early-types (E, SO, Sa) populate the upper right region, and 
late types (Sb, Sc, Irr) occupy the lower left. 

In order to test the completeness of this demarcation, 
we used the 62,190 g alaxies visually cla ssified as ellipticals 
in Galaxy Zoo data (|Lintott et al.l 2010). We cross-matched 
these to the control sample, selecting objects within 1" ra- 
dius. 34,606 matches were found, and 27,569 of these lie 
above tt*-r*>2.22 mag. Therefore, this demarcation gives 
~ 80% reliability (Fig(S] contours), suggesting ~ 20% of the 
early-type galaxy population may be excluded through ap- 
plying this criterion. 147,275 visually classified 'spiral' galax- 
ies from the Galaxy Zoo data are also present in our con- 
trol sample. 22% of these lie above •u*-r*>2.22, i.e. we ex- 
pect ~ 22% contamination in the early-type sample. Of the 
499,408 galaxies in the co ntrol sample, ~ 40% were classified 
as early-types, using the IStrateva et al.l colour selection, to 
leave 199,391 early-type galaxies in the control sample, mor- 
phologically selected in colour space at ~ 80% completeness 
and ~ 78% efficiency based on Galaxy Zoo visual analysis. 

We explore the 4,998 radio-loud AGN selected galaxies 
in colour space. Figure |S] shows the radio-loud AGN are pre- 
dominately early types, with > 70% lying above u*-r*>2.22, 
i.e. the radio-loud AGN sources are predominantly hosted 
by ellipticals as defined in colour space. The distribution in 
colour space is similar to that of the control sample, despite 
the AGN colour s possibly influencing th e light. In agreement 
with this result. [Griffith fc Sternl (|2010h found radio-selected 
AGN have a high incidence of being hosted by early-type 
galaxies. 

The 118 galaxies we flagged as star forming and re- 
moved from the radio-loud sample in £13.61 are shown as 

4 On a spectroscopically classified galaxy sample, early-types are 
recovered at >98% completeness and >83% reliability. 
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Figure 6. Colour-colour plot of radio-loud AGN (contoured, bin- 
size = 0.1 X 0.1 mag, 4 equally spaced levels, max density con- 
tour = 400/bin), and SF galaxies as defined in SI3.6I (crosses) The 
dashed line is u* — r* ^ 2.22 from IStrateva et al . (2001) and 
defines a morphological colour separator. 



crosses in Figure [6l and 98% of these sit below the line. 
We select the radio-loud AGN above the colour cut, to give 
3,516 'early-type' radio-loud AGN hosts. 



3.8 Redshift distributions 

At this stage we re-examined the redshift distributions of the 
samples, using the KS test. The probability of the control 
sample and the radio-loud AGN being drawn from the same 
parent sample in redshift was < 0.1%: the distributions in 



redshift differ (Figure 7(a) I. Evolutionary differences in host 
galaxy properties should be small over this entire redshift 
range, however we chose to re-sample the control sample 
redshift distribution to match that of the radio-loud AGN 
sample. 

We also examined the intrinsic r* band luminosity of 



the two samples. Figure 7(b) shows the radio sample is 
globally brighter in r* than the control sample. It has been 
well established that radio-loud AGN are hosted preferen- 
tially in the brightest and most massive elliptical galaxies 
jBest et al.ll2005bl ; iMauch fc Sadlerll2007l V Since we wanted 
to test whether radio-loud AGN harbour an excess of blue 
emission in their centres, we chose to re-sample the control 
sample to have the same r* band optical properties so as to 
avoid any possible correlation of colour gradient with galaxy 
luminosity. 

We generated a matched control sample by selecting the 
10 galaxies from the whole sample that lie closest to each 
radio-loud AGN in magnitude-redshift space. The final con- 
trol sample contains 35,160 galaxies, whose magnitude and 
redshift distributions are shown as dotted lines in Figure 
[7] This final selection should allow the distribution of ob- 
servables in the radio and control samples to be compared 
directly. 



3.9 Comparison radio sample (CR) 
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(b) Absolute r* magnitude distributions 

Figure 7. The fractional distributions of (a) redshift and (b) r* 
band luminosity. The solid black line shows the radio-loud AGN 
distributions. The red dashed line shows the control sample distri- 
bution prior to nearest neighbour matching in redshift-magnitude 
space. The blue dotted line shows the resultant sample randomly 
selected to match the redshift-magnitude distributions of the ra- 
dio sample. 



We imposed a redshift cut of 0.02 < z < 0.18 and re- 
stricted the redshift confidence parameter given by SDSS to 
> 95%, leaving 3,727 sources. We examined the CR sample 
radio powers based on NVSS and FIRST. Figure [5] shows 
the NVSS and FIRST luminosities for the 3,727 CR sources 
(grey crosses). Those with PvTghz > 10 23 WHz" 1 are high- 
lighted by blue squares (50%). We restricted the CR sam- 
ple to galaxies harbouring radio sources with P^XgHz > 
10 23 WHz -1 , which selects 1,847 (50%) of the sample. 

Of those 1,847 sources, 251 have FIRST derived powers 
< 10 23 WHz-\ and were therefore excluded in the radio 



The criteria applied to the 'radio' sample in £13.3113.81 were 
also applied to the 5,719 objects in the CR sample (see £|3.2[) . 



sample selection in £ j 3 . 4 1 

We removed 2 type 1 AGN and 48 star-forming galaxies 
confirmed at S/N > 3 from the CR sample, leaving 1,797 
radio-loud narrow-line AGN hosts. 4% of the sample are 
classified as type '1.9' AGN, with substaintial but not full 
obscuration of the central source. We applied the Strateva 
et al. colour cut (it*-r*>2.22 mag) detailed in ^3.71 to select 
1,295 early-type galaxies from the CR sample. 

The control sample was cross-matched with the FIRST- 
NVSS catalogue to within 2" and the 5 galaxies found were 
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Figure 8. FIRST vs NVSS derived luminosities for the 3,737 
sources in the candidate CR sample (grey crosses). Those with 
P i I 4GHz > 10 23 WHz" 1 are highlighted by blue squares (1,847). 
Of those, 251 sources with P^ 4 R J^ < 10 23 WHz -1 are shown as 
red crosses. These are the sources excluded in the radio sample's 
luminosity cut ( i]3.4|l . 

removed to ensure the comparison control sample is com- 
prised solely of 'radio-quiet' AGN elliptical hosts. 199,386 
galaxies remain in the CR-control sample. 

We selected the 10 galaxies from the control sample 
that lie closest to each CR radio-loud AGN in magnitude- 
redshift space. The final CR-control sample contains 12,950 
sources, matched to the CR sample in r*-band magnitude 
and redshift distributions. 

3.10 Summary of Final Samples 

Radio-loud 'early type' AGN (R-AGN) - 3,516 
sources 

After matching FIRST to the MGS (r* etro < 17.77 mag, 
0.02 < z < 0.18) within 2", we classified as radio-loud AGN 
hosts with Pi.4GHz > 10 23 WHz" 1 . Type 1 AGN exhibiting 
both broad Ha and H/3 emission lines were removed from 
the sample. Optical emission line ratios (if confirmed at 3cr) 
were used to remove star-fo rming galaxies using the demar- 
cation of lKauffmann et al.l (|2003f i. We then created a subset 
with u*-r* > 2.22 mag i.e. 'early-types' to remain consistent 
with the colour cut applied to the control sample. 
Control sample - 35,160 sources 

We removed all galaxies from SDSS-MGS (rp etro < 
17.77 mag, 0.02 < z < 0.18) with FIRST and NVSS coun- 
terparts. Type 1 AGN exhibiting both broad Ha and H/3 
emission lines were removed from the sample. Galaxies were 
constrained in colour to bias towards early-type morphology 
where a comparison with results of the Galaxy Zoo program 
suggests we have > 80% completeness and 78% efficiency. 
The resultant sample comprises radio-quiet, 'normal', pre- 
dominantly elliptical galaxies. We then selected a subsample 
to match the redshift and optical r*-band magnitude distri- 
butions of the R-AGN sample, thus removing any redshift 
bias and potential correlation of colour gradient with galaxy 
luminosity. 

Comparison radio (CR) - 1,295 sources 

We cross-matched the FIRST-NVSS catalogue derived by 



iKimball fc Ivezid l|2008l ) with the MGS (r* petro < 17.77 mag, 
0.02 < z < 0.18) to within 2", and classified radio-loud AGN 
hosts as those sources exhibiting P^Xghz > 10 23 WHz~ 1 . 
Type 1 AGN exhibiting both broad Ha and H/3 emission 
lines were removed from the sample. Optical emission line 
ratios (if confirmed at 3a) were us ed to remove star- f ormin g 
galaxies using the demarcation of lKauffmann et al.l l|200ot ). 
We then created a subset with u*-r* >2.22 mag i.e. 'early- 
types' to remain consistent with the colour cut applied to 
the control sample. 

CR-control sample - 12,950 sources 

All selection criteria as above for the control sample, we also 
removed 5 galaxies within the control sample that were also 
present in the CR sample. We then selected a subsample to 
match the redshift and optical r*-band magnitude distribu- 
tions of the CR sample. 

We therefore have two radio-loud AGN samples, both 
with a comparison control sample. The CR sample is derived 
using both radio catalogues, and is a brighter radio sam- 
ple (the NVSS flux density completeness limit is 2.5 mjy) 
for comparison with our deeper (> lmjy), larger but more 
incomplete R-AGN sample derived solely from the FIRST 
catalogue. 



4 ANALYSIS 

The surface brightness profile of an early-type galaxy can be 
described by the Sersic equation 

J(r) = I(r e )e- b ^ 1/n -^ (3) 

where r e is the effective radius (scale length) and I (re) is 
the corresponding effective surface brightness, b m 2n — 1/3 
is chosen so r e contains half the light in the galaxy, n w 4 for 
bright ellipticals d ecreasing to n ~ 2 as luminosity decreases 
l|Caon et al.lll993l ). SDSS chooses to fix n = 4, fitting the 
'de Vaucouleurs profile', truncating the profile beyond 7r e . 
The model fitted by SDSS has an arbitrary axis ratio and 
position angle and is convolved with a double- Gaussian rep- 
resentation of the PSF. The fitting then yields, among other 
properties; the effective radius, r e and error, Ar e . In general, 
the error for the u model is greater than other bands. It is 
noted that SDSS's fitting algorithm generates some weak 
discretization of model parameters, especially in r e (r) and 
r e (u). Objects with scale lengths lying in discrete bands in 
r e tend to have poorer goodness of fit. This is not included 
in the error estimates, Ar e , which were derived from count 
statistics on the image. In our work with the r e values, we 
consider Ar e rather than DeV_l (the likelihood of the model 
fit) due to our reluctance to remove data with poor de Vau- 
couleurs u band fitting. A poor goodness of fit to the de- 
Vaucouleurs profile may be indicative of a cuspy u band 
and removing such objects could bias results. 

SDSS denotes the f racDeV parameter as the fraction of 
the total flux contributed by the de Vaucouleurs component 
in a linear combination of a de Vaucouleurs and an exponen- 
tial model to find the best fit. It is noted that the likelihood 
values in the r*-band are intrinsically poor (all samples have 
mean values ~ 0.01), but the likelihood ratios still pick out 
reliable best-fit parameters. The fracDeV parameter (/ here- 
after) is correlated with the Sersic index; n = 1 corresponds 
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Figure 9. Normalised distributions of R for the radio-loud AGN 
(R-AGN) sample (solid) and the control sample (red dashed). 



to / = 0, n = 4 corresponds to / = 1 dKuehn fc Rvdeiill2005l . 
and references therein). Following I Vincent &: Rvden ( 20051 ). 
galaxies with 0.5 < / < 0.9 (2.0 <n< 3.3) are labeled 
"de/ex" galaxies and galaxies with / ^ 0.9 (n> 3.3) are la- 
beled "de" galaxies. We chose not to remove objects with 
/ < 0.5, but instead created subsets of each sample with 
/ > 0.5 for comparison. 

We investigate the colour structure of AGN host galax- 
ies principally through the distribution of R, which we define 
as R — r e (r)/r e (u), which is the ratio of r to u de Vau- 
couleurs effective radii. This should avoid any intrinsic scale 
differences in r and u between individual galaxies. 



4.1 Diversity of intrinsic distributions. 

The distributions of R for the R-AGN and control sample 
are shown in Figure[9] Both samples have a few (< 0.001%) 
strong outliers (R > 100) which significantly disturb the 
mean. The medians of the R-AGN and control distribu- 
tions are 0.837 and 0.793 respectively. Using a two-tailed 
Kolmogorov-Smirnov test we find the probability of the con- 
trol and the R-AGN sample being drawn from the same dis- 
tribution to be small (< 0.1%). 

Similarly, the CR sample's R distribution also shows a 
higher median compared to the CR-control sample (0.851 
and 0.797 respectively). The KS test also finds the proba- 
bility of the CR and CR-control sample being drawn from 
the same distribution in R is vanishingly small (< 0.1%). 



4.2 Maximum likelihood analysis 

The fractional errors in individual values of r e (it) are rel- 
atively large, which makes proving an intrinsic difference 
between the radio-loud AGN samples and their respective 
control sample distributions challenging. We used a maxi- 
mum likelihood analysis to compare the average values of R 
for the AGN and control sam ples, taking into a c count mea- 
surement errors according to iMaccacaro et all (|l988l ) and 
IWorralll l|l989h . in order to quantify the difference between 
the distributions of R. 

Under the assumption that both the AGN and control 
samples have normally-distributed intrinsic values of R, then 
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(b) CR and CR-control 

Figure 10. 90% (solid line) and 99% (dashed/dotted line) 
confidence contours of the mean ratio of scale lengths (R = 
r e (r)/r e (u)) and intrinsic dispersion for radio-loud AGN galaxies 
and 'normal' early types. The control and radio-loud 'early-type' 
AGN are distinct populations, as are the CR and CR-control sam- 
ples. 



the population mean, R, and intrinsic dispersion, a 
ther can be determined from individual values Ri ± a 
minimizing the function 
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Figure [10] shows the results of minimization of Eq [4] 
across R and a parameter space, for the two samples. The 
90% and 99% joint-confidence contours for R and a are given 
by S m in + AS for AS= 4.61 and 9.21 respectively. Table 
[2] shows the best fit R, a planes for both distributions. A 
clear separation between radio-loud 'early type' AGN and 
'normal' control galaxies can be seen for both the R-AGN 
and CR samples, and we conclude the two are not drawn 
from the same population at ^> 99% confidence. 

We applied the cut / > 0.5 in the r*-band, and re- 
examine the results for both radio-AGN samples. This cut 
will ensure the reliability of the de Vaucouleur's scale lengths 
used to derive R, and this subset is used as a comparison 
to the full samples to verify whether galaxies with shapes 
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Figure 11. 90% (solid line) and 99% (dashed/dotted line) 
confidence contours of the mean ratio of scale lengths (R = 
r e (r)/r e (u)) and intrinsic dispersion for the subset of galaxies 
with / > 0.5. The 1,851 radio-loud AGN galaxies and 17,338 
'normal' early types are shown in panel (a), the 575 CR galaxies 
and 6,364 CR-control galaxies are shown in panel (b). The con- 
trol and radio-loud 'early-type' AGN remain distinct populations. 
The smaller CR and CR-control sample lead to larger confidence 
contours that overlap, but these remain distinct at 89% confi- 
dence. 

that differ significantly from the de Vaucouleurs profile per- 
turb our main result. Table [2] also shows the fraction of each 
sample with 'good' de Vaucouleurs r*-band fits and their R 
values (see Figure . Although the values do not signifi- 
cantly differ from the full sample, the confidence contours 
now overlap in the smaller, but more complete, CR sample. 
However, the R-AGN and control samples are still seen to 
come from distinct populations. 



5 DISCUSSION 

Our results confirm that radio-loud AGN are hosted by 
brighter, bigger galaxies on average (see Figure^. Our tech- 
nique has shown the presence of an AGN is associated with 
an increase of the scale size of a galaxy in red light rel- 
ative to blue light. This suggests that radio galaxies have 
a bluer central bulge and/or more diffusely distributed red 



23.0 23.5 24.0 24.5 25.0 25.5 26.0 
Luminosity density (log,„(WHz~')) 

Figure 12. Radio luminosity density vs R for the CR sample. 
We found no correlation (> 99% certainty) between the two. 

light compared to their radio-quiet counterparts. With the 
analysis so far, we cannot tell whether this difference is from 
AGN driven or diffuse star-like light within the bulge, or 
more distributed red light. 

The CR sample is also significantly different in R rel- 
ative to the CR-control sample, so that our R-AGN sam- 
ple's lower completeness in extended sources does not signif- 
icantly affect the result. The CR sample dispersion is smaller 
than that of the R-AGN sample, whereas all four control 
samples have values of a that are similar. We applied a flux 
cut of > 4mJy to the R-AGN sample (61%), and found 
R = 0.80 ± 0.02 and a = 0.21 ± 0.01. The increased scatter 
about R is caused by fainter radio sources, which tend to be 
more heterogenous, and/or have underestimates of errors on 
r e values in more distant objects. 

The subset of CR galaxies constrained to have good r*- 
band de Vaucouleurs fits (/ > 0.5) does not contain enough 
radio galaxies to show strongly a significant separation from 
the CR-control sample (see Figure [TTT b)). We present dis- 
cussion predominantly on the R-AGN and control sample, 
which display a clear separation in both Figure UOf a) and 
lllf a). We concluded the two are not drawn from the same 
population at 2> 99% confidence. This sample is incomplete 
to radio sources with lobes > 5.4" and the population of 
galaxies with a weak radio AGN core but powerful lobes. 
However it is large due to the depth of the FIRST survey. 
The CR results for all Figures in the discussion are simi- 
lar to the results presented for the R-AGN sample. We find 
no significant deviation between the two in any subsequent 
par ameter examine d , exce pt in a, as discussed above. 

iMahabal et all l|1999l . hereafter MKM99) used a com- 
plementary ratid0 to demonstrate the ubiquity of inner blue 
components in a sample of only 30 radio galaxies and 30 con- 
trol galaxies, and argued that the blue light is due to star 
formation associat ed with the presence of a radio source (e.g. 
I Anton et al"1l2008l ). Figure 4 of MKM99 (power at 408 MHz 
vs log [r e (B)/r e (R)]) shows a trend for more powerful radio 
galaxies to have steeper colour gradients as indicated by the 
ratio of scale lengths. In contrast to this, when considering 
only the brightest (Pi.4GHz > 10 25 WHz- 1 ) radio sources m 

5 r e (B)/r e (R) 
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Sample 



R-AGN 
Control 
CR 



R 



/>0.5 



R 



f 



S 



0.76±0.02 
0.719±0.005 
0.78±0.02 



0.23±0.01 
0.206±0.003 
0.20±0.02 



CR-Control 0.731±0.007 0.204±0.006 



53% 
49% 
44% 
49% 



0.77±0.02 
0.722±0.006 
0.77±0.03 



0.25±0.02 
0.208±0.005 
0.19±0.03 



0.730±0.010 0.206±0.007 



Table 2. Best fit R and a corresponding to 90% error, for the R-AGN and CR sample, and their relative control samples. The fraction 
of each sample with r*-band / > 0.5 is shown, and the resultant mean and standard deviation, Rf and Of. 



the CR sample (for which the powers are more accurate) we 
find no correlation between radio power and R. We also find 
no correlation over the full range of radio powers and R as 
shown in Figure [121 the Spearman's rank correlation result 
is ±0.07 or less between Pi. 4 ghz and R. We infer that a high 
value of R (and possible blue central bulge) is not associated 
with the power of the radio source. 

There are several possible scenarios to be considered to 
explain the larger R in the R-AGN sample. 

5.1 AGN light 

Central point-like it* light may be indicative of a blue AGN. 
We investigated the nature of the R-AGN excess R, using the 
goodness of the deVaucouleurs fit. The SDSS deVaucoulcurs 
fitting procedure returns the likelihood (between and 1) 
associated with the model from the \ 2 n t. If the u*-band 
goodness of fit for the R-AGN sample approaches 1 for low 
R and for high R then this is indicative of an AGN/core 
component (as the AGN point-like component would per- 
turb a deVaucouleurs fit) driving a blue excess and high R. 
If no structure is seen in R vs deVaucouleurs u*-band good- 
ness of fit in either sample, then any blue excess is more 
likely driven by some sort of diffuse starlight in the central 
bulge. 

The w*-band data for both the R-AGN and control sam- 
ples are generally well fit by a deVaucouleurs profile (the 
R- AGN/control median likelihood values are 0.86 and 0.83 
respectively), and there is no trend for the u*-band or r*- 
band fits to be worse at large R in the R-AGN than the 
control sample. We can conclude that the high R in the R- 
AGN sample is not due to any point-like AGN component. 

5.2 Radio AGN triggering star formation 

It is generally agreed that AGN feedback regulates the sup- 
ply of cold gas for star formation by supplying energy to the 
ISM and IGM. Some radio-loud star-forming AGN show a 
connection between the suppression of star formation and 
the strength of th e radio jets, by heatin g and expelling the 
surrou nding gas l|Nesvadba et alj 120081 ) . ISchawinski et afl 
(2009) looked at a sample of low redshift SDSS early- 
type galaxies for which late-time star formation is being 
quenched. They found that molecular gas disappears less 
than 100 Myr after the onset of accretion onto the central 
black hole. These galaxies were not associated with radio 
jets, but show that low-luminosity AGN episodes are suffi- 
cient to suppress residual star-formation in early-type galax- 
ies. 

The residual fractional number density distribution of 
scale lengths are shown in Figure [131 for the control sample 
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Figure 13. Residual fractional number density distributions (R- 
AGN - control), in 2 kpc bins, for the u*-band and r*-band scale 
radii. There is a higher fraction of blue cores in the R-AGN sample 
as compared to the control sample at r e ~ 10 kpc. This plot does 
not demonstrate the correlation between the red and blue scale 
radii (see Figure HH)! . The measurement errors are larger in the u* 
band which may account for the higher variability in the residuals. 

subtracted from the R-AGN sample. There is a clear ex- 
cess of blue cores (0.1 — 10 kpc) in the R-AGN sample as 
compared to the control sample. There are also more small 
r e (r*) in the R-AGN sample, ranging from 4 — 10 kpc. This 
plot does not demonstrate the correlation between the red 
and blue emission, which is shown in Figure fT4l but it does 
demonstrate that the presence of a radio AGN seems not to 
have suppressed star-formation in the central regions of its 
host galaxy. 

Figure [141 plots r e (u*) against r e (r*) for the R-AGN 
sample (black) and the control sample (red dashed). The 
dotted - dashed line denotes the locus r e (r*) = 0.719r e (w*) 
as defined by the control sample (see Table [2j . 

At r e (r*) > 0.719r e («*), the surface brightness in u* 
increases more rapidly toward the centre of a galaxy than 
in r* as compared to the mean of the control sample, i.e. 
half the blue light from the galaxy is contained in a smaller 
region than half the red light, implying the galaxy becomes 
bluer inward. 

The distribution of number densities in Figure [T4l shows 
radio galaxies appear to become bluer inwards more often 
than the control galaxies (66% of the R-AGN sample are 
above the locus r e (r) = 0.719 r e (u), compared to 60% of the 
control sample). 

There are fewer blue smaller cores in the control sam- 
ple as compared with the R-AGN sample (see also Figure 
I13[) , implying there may be a blue excess near the AGN 
core for some of the radio galaxies. However, the higher R 
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Figure 14. Normalised number densities of r e (u*) and r e (r*) for 
the R-AGN (solid) and control sample (red dashed). The dotted - 
dashed line is the locus r e (r*) = 0.719r e («*), galaxies above this 
line become bluer inward. Contours show the fractional number 
densities of each sample with levels at 0.004, 0.008 and 0.011. 
There is a larger fraction of the R-AGN sample above the line, 
possibly denoting a tendency to become bluer inward more often 
than the control sample. 



values of the R-AGN sample also seems to be driven by a 
higher scale length in r* , seen in the difference between the 
inner two contour levels of each sample (see also i|5.4[) . so 
we cannot definitively conclude that the increased R in the 
R-AGN sample is due to star formation within the central 
few kpc. Star-formation within the central few kpc would 
contradict feedback models which predict the suppression 
of star formation near an AGN. 



5.3 Star formation in the bulge 

IShabala et alj (|2008l ) found that radio sources in massive 
hosts are re-triggered more frequently than their less massive 
counterparts, suggesting that the onset of an AGN quiescent 
phase is due to fuel depletion. AGN activity is therefore 
promoted by an increase in gas in the centres of the galaxies, 
which may imply a link between the AGN radio phase and 
star formation in the bulge (< lOkpc) of the host galaxy. 

It has been suggested that AGN activity and a major 
episode of star formation in radio-loud galaxies is triggered 
by the accretion of gas during major mergers and/or tidal 
interactions. However AGN activity is initiated later in the 
merger event than t h e starburst fe.g.lSchawinski et al.ll2007l ; 
lEmonts etafl l200rj ; iTadhunter et all 120051 ). Our spectro- 
scopic selection rules out major merger events (see H2.1.1I) . 
but residual star formation may still be present on global 
scales. 

Observation ally, there is a s trong link between AGN 
and starbursts ( Shin et al.l 1201 ll . and references therein). 
iKauffmann et all found powerful optical AGN (as 

classified by the strength of the [OIII] emission L[om] > 
10 7 Lq) predominantly reside in 'young bulges'. Recent star 
formation can provide up to 25 — 40% of the optical/UV 
conti nuum in radio galaxies at low and inte rmediate z e.g. 
(e.g.. lHolt et al.ll2007j ; ITadhunter et alj|2005h . 

We found a higher fraction of R-AGN galaxies have 
r e (u*) < 10 kpc as compared to the control sample (see Fig- 



ure !13[) . perhaps indicative of star formation in the outskirts 
of the bulge, fuelled by gas expelled from the central regions 
by the AGN. 



5.4 Diffuse red emission 

An alternative interpretation of a larger R is that the R- 
AGN sample has more distributed red light. Radio- loud 
'early-type' galaxies predominantly reside in elliptical galax- 
ies, which are well k nown to be redder than late- types (e.g., 
IStrateva et al.l [200 lh . Diffuse red emission would cause a 
larger deVaucouleurs scale radius in the r* band and an 
increased R. 

Figure [T4l suggests the increased R in the R-AGN sam- 
ple may also be attributed to more diffuse red light; between 
r e (u*) ~ 4— 10 kpc, r e (r*) is on average higher for the radio 
population. This is also seen in the lower panel of Figure fTBI 

Our results show a difference in the distributions of red 
and blue light in R-AGN and normal 'early-type' galaxy 
populations, though a high value of R is not a property 
of every active galaxy. We found the higher R in the R- 
AGN sample is contributed to by a higher fraction of radio 
galaxies harbouring small blue cores but also an increase in 
the numbers of radio galaxies with more diffuse red light as 
compared to the control sample. Given that star formation 
proceeds over a longer timescale than radio activity, this 
disfavours the idea that all galaxies undergo short bursts 
of radio activity, but rather implies that a subset have the 
predisposition to become radio-loud. 



6 CONCLUSION 

We cross-matched low-redshift (0.02 < z < 0.18) data from 
the SDSS MGS (r* etro < 17.77 mag) and FIRST to within 
2", deriving a radio sample of galaxies at > 99% efficiency 
and > 72% completeness. Radio luminosities were in the 
range 10 23 - 10 25 WHz" 1 . Type 1 AGN were removed from 
the sample (identified via Ha and H/3 broad line characteris- 
tics), 4% of the sample are classified as type '1.9' AGN, with 
substantial but not full obscuration of the central source. 
Contamination from SFGs (identified via optical emission 
line ratios) in the radio sample is expected at ~ 11%. A con- 
trol sample was defined from SDSS sources with no match 
to a FIRST source within 2" of their optical core, providing 
a sample of > 99% efficiency and > 92% completeness. At 
80% reliability, the demarcation u* — r* > 2.22 mag selected 
'early-type' galaxies in both samples. Samples were matched 
in r*-band magnitude and redshift distributions and final 
sample sizes were 3,516 radio-loud AGN galaxies (R-AGN) 
and 35,160 control galaxies. 

We also created a complementary flux-limited sample 
through cross-matching with NVSS (the CR sample). The 
same cuts were applied to derive a radio-loud 'early-type' 
AGN sample, except we used the more all-encompassing 
NVSS flux estimates to cut in radio luminosity, thereby re- 
taining the populations of galaxies with weak/no core ra- 
dio emission but bright, extended radio lobes. This sam- 
ple had higher completeness for comparison with the R- 
AGN sample. A control sample was defined from the SDSS 
sources with no match to a NVSS-FIRST source within 2" of 
their optical cores. We further considered only sources where 
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fracDeV (f) > 0.5 to restrict all four samples to having good 
r*-band deVaucouleurs fits. 

We investigated the colour structure of AGN host galax- 
ies through R, the ratio of r* to u* de Vaucouleurs effective 
radii and used maximum likelihood analysis to quantify the 
degree of difference in the distribution of R between sam- 
ples. We concluded the radio (R-AGN/CR) samples are not 
drawn from the same population as their radio-faint control 
samples at S> 99% confidence: the presence of an AGN in- 
creases the scale size of a galaxy in red light relative to blue 
light, on average. 

Our result does not appear to be driven by the presence 
of blue AGN in the radio-loud samples since the goodness 
of the de Vaucouleurs fits does not become worse as R in- 
creases. We found no structure in R vs u*-band goodness of 
fit in either radio sample. We found an excess of blue cores 
in radio-loud galaxies as compared to radio-quiet, 'early- 
type' galaxies, implying the increased R may be due to star 
formation in the central few kpc, in contrast with feedback 
models which predict the suppression of star formation near 
an AGN. Spectroscopic selection of our samples rules out 
major merger events, as starbursts can be triggered by the 
accretion of gas/tidal interactions. We also found radio AGN 
hosts to have larger red scale lengths in relation to their blue 
light and note this to be a contributing factor in an increased 
R. We cannot definitively discern whether a small blue core 
or larger distribution of red light is the driving factor in this 
result. 

Given the longer timescale for star formation than ra- 
dio activity, our results imply a subset of galaxies have the 
predisposition to become radio-loud, rather than all galax- 
ies undergoing bursts of radio activity at some stage in their 
lifetimes. 



ACKNOWLEDGEMENTS 

EM gratefully acknowledges support from the UK Science 
and Technology Facilities Council and thanks Luke Davies 
and James Price for helpful comments. We thank the anony- 
mous referee for helpful and insightful comments that have 
resulted in an improved analysis. In undertaking this re- 
search, we m ade extensive use of the TOPCAT software 
i|Tavloij|2005T ). Funding for the SDSS has been provided 
by the Alfred P. Sloan Foundation, the Participating In- 
stitutions, the National Science Foundation, the U.S. De- 
partment of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Coun- 
cil for England. This research makes use of the NVSS and 
FIRST radio surveys, carried out using the National Ra- 
dio Astronomy Observatory Very Large Array: NRAO is 
operated by Associated Universities Inc., under cooperative 
agreement with the National Science Foundation. 



REFERENCES 

Abazajian K. N., Adelman-McCarthy J. K., Agiieros M. A., 
Allam S. S., Allende Prieto C, An D, Anderson K. S. J., 
Anderson S. F., Annis J., Bahcall N. A., et al. 2009, ApJS, 
182, 543 



Anton S., Browne I. W. A., Marcha M. J., 2008, AAP, 490, 
583 

Antonucci R., 1993, ARAA, 31, 473 

Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP, 
93, 5 

Becker R. H., White R. L., Helfand D. J., 1995, ApJ, 450, 
559 

Best P. N., Kauffmann G., Heckman T. M., Brinchmann 
J., Chariot S., Ivezic Z., White S. D. M., 2005, MNRAS, 
362, 25 

Best P. N., Kauffmann G., Heckman T. M., Ivezic Z., 2005, 

MNRAS, 362, 9 
Blanton M. R., Dalcanton J., Eisenstein D, Yanny B., Ya- 

suda N., York D. G., 2001, AJ, 121, 2358 
Blanton M. R., Roweis S., 2007, AJ, 133, 734 
Caon N., Capaccioli M., D'Onofrio M., 1993, MNRAS, 265, 

1013 

Condon J. J., 1989, ApJ, 338, 13 

Condon J. J., Cotton W. D, Greisen E. W., Yin Q. F., 
Perley R. A., Taylor G. B., Broderick J. J., 1998, AJ, 115, 
1693 

Emonts B. H. C, Morganti R., Tadhunter C. N., Holt J., 
Oosterloo T. A., van der Hulst J. M., Wills K. A., 2006, 
aap, 454, 125 

Fanaroff B. L., Riley J. M., 1974, MNRAS, 167, 31P 
Fukugita M., Ichikawa T., Gunn J. E., Doi M., Shimasaku 

K., Schneider D. P., 1996, ApJ, 111, 1748 
Gonzalez Delgado R., Tadhunter C, Perez E., 2006, As- 

tronomische Nachrichten, 327, 159 
Gonzalez-Perez V., Castander F. J., Kauffmann G., 2010, 

MNRAS 

Graham A. W., Driver S. P., Petrosian V., Conselice C. J., 
Bershady M. A., Crawford S. M., Goto T., 2005, AJ, 130, 
1535 

Graham A. W., Erwin P., Trujillo I., Asensio Ramos A., 

2003, AJ, 125, 2951 
Griffith R. L., Stern D, 2010, AJ, 140, 533 
Gunn J. E., Hogg D, Finkbeiner D, Schlegel D, 2001 
Holt J., Tadhunter C. N., Gonzalez Delgado R. M., Inskip 

K. J., Rodriguez Zaurin J., Emonts B. H. C, Morganti 

R., Wills K. A., 2007, MNRAS, 381, 611 
Hutchings J. B., Neff S. G., Weadock J., Roberts L., Ryn- 

eveld S., Gower A. C, 1994, AJ, 107, 471 
Ivezic Z., Menou K., Knapp G. R., Strauss M. A., Lupton 

R. H., Vanden Berk D. E., Richards G. T., Tremonti C, 

Voges W., White R. L., Yanny B., York D. G., 2002, AJ, 

124, 2364 

Kauffmann G., Heckman T. M., Tremonti C, Brinchmann 
J., Chariot S., White S. D. M., Ridgway S. E., Brinkmann 
J., Fukugita M., Hall P. B., Ivezic Z., Richards G. T., 
Schneider D. P., 2003, MNRAS, 346, 1055 

Kimball A. E., Ivezic Z., 2008, ApJ, 136, 684 

Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, 
C. A., & Trevena, J. 2001, ApJ, 556, 121 

Kuehn, F., & Ryden, B. S. 2005, ApJ, 634, 1032 

Lin Y., Shen Y., Strauss M. A., Richards G. T., Lunnan 
R., 2010, APJ, 723, 1119 

Lintott C, Schawinski K., Bamford S., Slosar A., Land K., 
Thomas D., Edmondson E., Masters K., Nichol R., Rad- 
dick J., Szalay A., Andreescu D., Murray P., Vandenberg 
J., 2010, ArXiv e-prints 

Lu Y., Wang T., Zhou H., Wu J., 2007, AJ, 133, 1615 



14 E. J. A. Mannering et al. 



Maccacaro T., Gioia I. M., Wolter A., Zamorani G., Stocke 

J. T., 1988, ApJ, 326, 680 
Mahabal A., Kembhavi A., McCarthy P. J., 1999, ApJL, 

516, L61 

Masci F. J., Cutri R. M., Francis P. J., Nelson B. O., 
Huchra J. P., Heath Jones D., Colless M., Saunders W., 
2010, PASA, 27, 302 

Mauch T., Sadler E. M., 2007, MNRAS, 375, 931 

Moric I., Smolcic V., Kimball A., Riechers D. A., Ivezic Z., 
Scoville N., 2010, APJ, 724, 779 

Nesvadba N. P. H., Lehnert M. D., De Breuck C., Gilbert 
A. M., van Breugel W., 2008, AAP, 491, 407 

Oke J. B., Gunn J. E., 1983, ApJ, 266, 713 

Osterbrock D. E., 1989, Astrophysics of gaseous nebulae 
and active galactic nuclei 

Sadler E. M., Jackson C. A., Cannon R. D., Mclntyre V. J., 
Murphy T., Bland-Hawthorn J., Bridges T., Cole S., Col- 
less M., Collins C, Couch W., Dalton G., 2002, VizieR 
Online Data Catalog, 732, 90227 

Schawinski K., Lintott C. J., Thomas D., Kaviraj S., Viti 
S., Silk J., Maraston C, Sarzi M., Yi S. K., Joo S., Daddi 
E., Bayet E., Bell T., Zuntz J., 2009, ApJ, 690, 1672 

Schawinski K., Thomas D., Sarzi M., Maraston C, Kaviraj 
S., Joo S.-J., Yi S. K., Silk J., 2007, MNRAS, 382, 1415 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Shabala S. S., Ash S., Alexander P., Riley J. M., 2008, 

MNRAS, 388, 625 
Shin M., Strauss M. A., Tojeiro R., 2011, MNRAS, 410, 

1583 

Stoughton C, Lupton R. H. e. a., 2002, AJ, 123, 485 
Strateva I., Ivezic Z., Knapp G. e. a., 2001, AJ, 122, 1861 
Strauss M. A., Weinberg D., Lupton R., et al. 2002, ApJ, 
124, 1810 

Tadhunter C, Robinson T. G., Gonzalez Delgado R. M., 
Wills K., Morganti R., 2005, MNRAS, 356, 480 

Taylor, M. B. 2005, Astronomical Data Analysis Software 
and Systems XIV, 347, 29 

van Breugel W., Fragile C, Croft S., de Vries W., Anninos 
P., Murray S., 2004, in Storchi-Bcrgmann T., Ho L. C, 
Schmitt H. R., eds, The Interplay Among Black Holes, 
Stars and ISM in Galactic Nuclei Vol. 222 of IAU Sym- 
posium, Jet-induced star formation: Good news from big, 
bad black holes, pp 485-488 

Vincent, R. A., & Ryden, B. S. 2005, ApJ, 623, 137 

Worrall D. M., 1989, in J. Hunt & B. Battrick ed., Two 
Topics in X-Ray Astronomy, Volume 1: X Ray Binaries. 
Volume 2: AGN and the X Ray Background Vol. 296 of 
ESA Special Publication, Population properties of AGN 
in the soft X-ray band, pp 719-726 

Yan R., Ho L. C, Newman J. A., Coil A. L., Willmer 
C. N. A., Laird E. S., Gcorgakakis A., Aird J., Barmby P., 
Bundy K., Cooper M. C, Davis M., 2010, ArXiv e-prints 

Yan R., Newman J. A., Faber S. M., Konidaris N., Koo D., 
Davis M., 2006, ApJ, 648, 281 

York D. G., Adelman J., Anderson Jr. J. E., Anderson 
S. F., Annis J., Bahcall N. A., Bakken J. A., Barkhouser 
R., Bastian S., et al 2000, AJ, 120, 1579 

Yun M. S., Reddy N. A., Condon J. J., 2001, ApJ, 554, 803 



